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Parameterized frameworks for modified gravity are potentially useful tools for model-independent 
tests of General Relativity on cosmological scales. The toy model of an Einstein-de Sitter (EdS) 
universe provides a safe testbed in which to improve our understanding of their behaviour. We 
implement a mathematically consistent parameterization at the level of the field equations, and use 
this to calculate the evolution of perturbations in an EdS scenario. Our parameterization explicitly 
allows for new scalar degrees of freedom, and we compare this to theories in which the only degrees 
of freedom come from the metric and ordinary matter. The impact on the Integrated Sachs- Wolfe 
effect and canonically-conserved superhorizon perturbations is considered. 



I. INTRODUCTION 

A successful theory of modified gravity has proved hard 
to find. Most theories begin by postulating new physi- 
cal principles, for example, the existence of fundamental 
fields with certain symmetry properties [T] or additional 
dimensions [21 E]; see [H [S] for comprehensive surveys 
of the current literature. By introducing such principles 
on which to base a theory we are immediately select- 
ing specific directions in theory space to investigate. An 
alternative approach is to cautiously explore outwards 
from the corner of theory space that we understand best, 
that is, General Relativity (GR). One way to implement 
this strategy is to construct a parameterized framework 
that allows for small deviations from GR in a cosmolog- 
ical context, akin to the Parameterized Post-Newtonian 
framework that has been used to test GR extensively 
within the Solar System 

In order to extract maximum benefit from this ap- 
proach we need to develop a sense of how the parameters 
we use impact the growth of structure. Historically, our 
theories of structure formation have been developed using 
matter-dominated cosmological models [TUHT^ . Though 
not appropriate for realistic calculations, a toy Einstein- 
de Sitter (EdS) cosmology remains an immensely useful 
testing ground for new theories; here, we can gain a qual- 
itative understanding of perturbation evolution whilst 
our knowledge of the background expansion remains on 
a firm footing. Therefore, maintaining the principle of 
caution advertised above, we will consider the effects of 
parameterized gravity on an EdS universe (which con- 
tains pressureless matter only) . This will be an excellent 
approximation to the matter epoch of the real universe, 
and has the added advantage that analytic solutions are 
achievable in some cases due to the simple properties of 
pressureless matter. 

Despite the wide variety of modified gravity theories 
present in the literature, a survey of field equations re- 
veals some common features. In many theories the grav- 
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itational constant appearing in the Poisson equation ac- 
quires a time-dependence and/or scale-dependence |15| . 
and the Newtonian gravitational potential 'I' and curva- 
ture perturbation $ (in the conformal Newtonian gauge) 
are not equal as they are in GR. Parameterized frame- 
works for modified gravity are usually constructed to in- 
corporate such properties. We can ask what typical ef- 
fects these generic features might have on observables: 
for example, in addition to affecting the growth rate of 
structure, distortion of gravitational potentials will im- 
print secondary anisotropics on the CMB. However, we 
will assume that modifications to GR must be negligi- 
ble at very early times, to avoid significantly impacting 
the primary GMB and the sensitive reaction rates of Big 
Bang Nucleosynthesis J,6 18J; see [12] for consideration 
of scalar-field models that are non-negligible at the time 
of recombination. 

This paper investigates the evolution of cosmologi- 
cal perturbations in the parameterized framework imple- 
mented in [20^23) . for an EdS background. |jll] intro- 
duces the necessary formalism in the context of theories 
that are constructed purely from metric quantities. In 
pil| we use this framework to calculate the evolution of 
density perturbations on intermediate and large scales. 
Two other quantities of interest are also calculated: the 
growth function f{z), and the Integrated Sachs- Wolfe ef- 
fect that is induced in such gravity theories (note that 
this is not the same as the late-time Integrated Sachs- 
Wolfe effect that occurs in a A-dominated era). 

In §IV| we take the first steps towards implementing 
a similar treatment of gravitational theories that intro- 
duce additional scalar degrees of freedom. New degrees of 
freedom cause a significant increase in complexity which 
renders a general, model-independent calculation almost 
impossible, at least analytically. We will implement a 
parameterized effective fluid approach known as 'Gener- 
afized Dark Matter' (GDM) [Stj to facilitate the treat- 
ment of these additional scalars, and consider the case 
in which the effective fluid has a negligible equation of 
state. This approximation excludes some classes of the- 
ories from our analysis, such as f{R) gravity, but is ap- 
plicable in other cases; we drop the restriction again in 
fV] . In ^[V] we consider a perturbation that is conserved 
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on super-horizon scales in GR, and ask under what con- 
ditions this fact remains true in parameterized modified 
gravity. The conclusions of this work are presented in 

m 

The purpose of this paper is to develop an understand- 
ing of parameterized gravity, not to pursue accurate cal- 
culations for the real universe. Hence the plots and trial 
parameter values used here are intended to be illustrative 
rather than realistic. 



II. PARAMETERIZATION OF METRIC-ONLY 
GRAVITY THEORIES 

No standard parameterization of modified gravity cur- 
rently exists. A common choice is to introduce a free 
function that describes any time- or scale- variation of the 
gravitational constant in the Poisson equation [251 - 130] : 



(1) 



where k = SttGq, Gq is the canonical value of New- 
ton's gravitational constant, A = 5 -I- 3719 is a comov- 
ing gauge-invariant density perturbation and 6 is the ve- 
locity potential of a fluid defined by Vi — ViO. Often a 
second free function is used to describe the ratio of the 
two conformal Newtonian potentials, ri{a,k) = ^/vP. By 
redefining ( = (1 — 1/?^) this can be rewritten: 



$ - * = C(a,/fc)«' 



(2) 



The arguments of the 'modified gravity functions' 
(MGFs) /i and ( will be suppressed hereafter. We will 
call this form of the slip relation parameterization B. In 
|22j it was argued that a parameterization of this type is 
only applicable in the quasistatic regime, and that in the 
case of purely metric theories it implicitly corresponds 
to higher-derivative theories. This is true if one assumes 
that eqn.Q is an exact 'template' for the slip relation 
of a modified gravity theory, which means that ^ can 
only be a function of homogeneous background quanti- 
ties. If instead one is prepared to let C be a function 
of non-homogeneous environmental variables and initial 
conditions then the slip relation may not necessarily im- 
ply a higher-derivative theory; it becomes impossible to 
ascertain the derivative order of the theories being pa- 
rameterized without further information [5T| . 

In [35] an alternative format was suggested, in which 
the order of derivatives in the field equations is made 
explicit. In this alternative parameterization a metric 
theory with a ACDM background incurs extra constraint 
equations (see Table IT] in the appendix) , which force the 
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slip relation to be |20ri21) : 

<^> - -qj = C'^ + il^—^^ 
rtp 



MGF fi appears in both the Poisson and slip equations, 
unlike parameterization B. This apparently small differ- 
ence leads to some degree of ambiguity in the interpreta- 
tion of current constraints on the MGFs [23] . 

However, the purpose of this paper is not to discuss 
the subtleties of parameterization choice at length; we 
wish to keep our results as general as possible. We can 
treat both of the parameterizations simultaneously by 
adopting the following Poisson and slip equations: 



^J■s-^ 



$ - * = 



'H^J.s 



(4) 
(5) 



This slip relation will be termed parameterization A. The 
key feature to note here is that in parameterization A the 



To recover parameterization B we set /i^ = 1 but keep pp 
general. To recover parameterization A we set = pp. 
Note that ps is related to the function g in [22 by 
Ps = (1 — g)^^. Throughout this section we will some- 
times leave results expressed in terms of the three func- 
tions pp, Ps and C- We wish to emphasize from the out- 
set that there are only really two independent functions, 
and all expressions should be evaluated in either the A- 
type or B-type instance. We write our expressions in this 
general format because can be instructive to see whether 
the modified terms have their origin in the Poisson equa- 
tion (indicated by the presence of pp) or the slip relation 
(indicated by Ps and C)- 

What theories map onto equations Q and ([5|? The 
answer is 'very few', which is a cause for concern given 
that the above forms are often used to obtain constraints 
on modified gravity from current data. To map exactly 
onto these parameterizations a theory must stem from 
an action that is constructed only from curvature invari- 
ants and leads to fields equations that contain at most 
second- or third-order time derivatives (for parameteri- 
zations A and B respectively). By the immense power of 
Lovelock's theorem [32l [33] such a theory can only differ 
from GR if it introduces either nonlocality or spacetimes 
of dimension greater than four. This is a very restricted 
class of theories, though examples do exist [5iH5S] . 

The limitations described above can be relaxed if one 
is satisfied with an approximate correspondence between 
theories and parameterization, rather than an exact one. 
For example, it is frequently assumed that the perturbed 
Einstein equations will retain the form of eqns.Q and 
([2| in theories with additional scalar degrees of freedom. 
This cannot be exactly true; any new scalar coupled to 
gravity will modify the zeroth-order Einstein equations 
in some way, and we expect to see perturbations of the 
new field appearing in the linearized Einstein equations. 
However, the form of eqns.Q and ^ is retained within a 
limited range of distance scales for some theories [12 1571 - 

Soj. 

To avoid such approximations we will proceed by tak- 
ing eqns.(|4]) and ([s]) at 'face value', i.e. assuming that 
there are no new scalar degrees of freedom hidden be- 
hind them. This is the assumption that is implicitly be- 
ing made if equations such as ([T]) and 
in an Einstein-Boltzmann solver 



21 are implemented 
and used to 
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generate ISW and matter power spectra. In { IV we will 
introduce an extended parameterization that attempts to 
account for the additional scalars explicitly. 



III. DENSITY PERTURBATIONS 

In this section we will consider how the growth of cold 
dark matter (CDM) density perturbations in an EdS uni- 
verse is influenced by the MGFs. This is a model for the 
matter-dominated epoch of the real universe. The growth 
of structure during an epoch in which dark energy is also 
relevant was investigated in |21j . using a specific ansatz 
for and 

We will assume that any time-variation of the MGFs 
during the matter-dominated epoch must be very small 
in order to prevent them from evolving to a region of 
parameter space that would cause conflict with observa- 
tions. This is not guaranteed, for example, in [23] we 
found that cancellations between MGFs can lead to ob- 
servables very close to the predictions of ACDM. Hence 
even very radical departures from GR can be accommo- 
dated by current data in finely-tuned situations; how- 
ever, for the purposes of this paper we will assume that 
our present universe does not correspond to such a case. 
We will therefore take the time derivatives of /xp, /u^ and 
( to be negligible in comparison to the rate of evolution 
of other variables. We intend to relax this restriction in 
future work. 

Although we are neglecting their time-dependence, /ip 
and C may still contain scale-dependence. They are di- 
mensionless functions, but scale-dependence can appear 
as a ratio to some special scale that arises in a given the- 
ory, i.e. k/k^,. An example of such a privileged scale is 
the Compton wavelength of the scalaron in f{R) gravity 

The fluid conservation equations for CDM energy den- 
sity and momentum perturbations are given by: 



5m 



'k^OM 



3$ 



(6) 
(7) 



where 9m is the velocity potential. Differentiating 
eqn.([6]) and combining with cqn.([7]) leads to a second- 
order equation for 5m'- 
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1-L5m - 3$ - + fc2* = 



(8) 



This is the same as found in GR. However, differences 
from GR arise when we use a non-trivial slip relation to 
eliminate . Substituting eqn.([5]) into eqn.Q: 



3 a/- 



•H(5m-3$-3-H$ 



k"^ f fi 
3^2 



1 



+k^{i-0^ = 

(9) 



Before studying the behaviour of this equation it is useful 
to delineate a hierarchy of distance scales: 



1. The nonlinear scale on which clusters and galaxies 
form. 

2. The quasistatic scale on which time derivatives of 
perturbations can be neglected in comparison to 
their spatial derivatives. 

3. Larger scales on which the above approximation is 
no longer valid, but are still well within the horizon. 

4. Scales that are greater than our observable horizon. 

We will consider the solutions of eqn.(|9| in regions 3 and 
4. In fact it is possible to derive a single equation for 
$ that is valid in both regions 3 and 4, and then use 
the Poisson equation to relate its solutions to 5m (we 
thank C. Skordis for pointing this out). We will use an 
equivalent method that is simpler but a little less elegant. 



A. Subhorizon scales 

In region 3 described above we can approximate 
'H/k ^ 1 and Am ~ 5m since \9m\ ^ \'^M\/\k\ is small. 
We use derivatives of eqn.Q to eliminate $ from eqn.(|9|: 



:>M 



nsM 



2 fis 



(10) 



M 



1-C 



(m. - 1) 



= 



where we have used the result H = —\T-P in an EdS 
universe. Writing the solutions of this equation in the 
form 



(11) 



where and N are constants and the power-law in- 
dices are: 



3^j 



(m. - 1) 



(12) 



±- 



2 

9^(m. - If - ii)— + 6Aip(9 -40 + 1 



In parameterization A this reduces to 



"a = 1 - oA* 



± ^VV + 12/i(3-2C)-20 



2^y 2 

whereas in parameterization B it becomes 
4 = -^±^v/l + 24Atp(l-C) 



(13) 



(14) 



It can be verified that in the limit /is = /ip = 1 and 
C — eqn.(12) recovers the GR result 5m oc a. We can 



see immediately that in both parameterizations there is 
a term — 24/ipC that leads to degeneracy between the ef- 
fects of the individual MGFs. Note that C, only appears 
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within this degenerate combination, so it cannot signif- 
icantly impact growth if fj,p is small. We note that if a 
/i-like MGF is implemented in the Poisson equation for 
^' instead of $ then this degeneracy does not arise in the 
parameterization B case ^28, ,39^ i42j . 



When 



are imaginary the solutions for Sm are 



damped oscillations. However, since our calculation has 
neglected the effects of baryons or radiation this oscil- 
latory behaviour simply indicates unphysical solutions 
rather than anything meaningful. To have at least one 
growing mode we need n'^ to be positive, for which the 
relevant condition is: 



M.(2-C) > 1 



(15) 



We expect that an approximate version of this bound 
should be obeyed in the real universe, in order to repro- 
duce the observed matter power spectrum - see ^ III C 



However, in the real universe the hard bound of eqn.( 15 ) 



will be softened by contributions to growth from the ra- 
diation and A-dominated eras. Note that there is no re- 
striction that prevents ( from adopting negative values. 

We wish to understand the physical mechanisms 
through which the MGFs are exerting their influence on 
small scales. We can get a feel for this by thinking about 
eqn.(ll) in the context of a simple mechanical system. 



The last term on the righthand side represents a time- 
dependent forcing that drives the collapse of density per- 
turbations. The 6m term is analogous to a frictional 
force, which in familiar physical situations always acts to 
oppose motion; its magnitude decreases with time due to 
the factor of Ji. Since we are parameterizing around a 
ACDM background the evolution of H is unaffected by 
the MGFs, and hence cannot be contributing to devia- 
tions from GR. 

The overall magnitude of the driving term is controlled 
by Hp. This intuitively makes sense - if we increase 
the gravitational coupling strength then we expect struc- 
tures to collapse faster. Less intuitive is the appearance 
of lis and C in the driving term, which have the abil- 
ity to change its sign. We will assume throughout that 
MPj Ms > always to maintain agreement with our phys- 
ical notion of attractive gravity, but note that there is no 
such restriction on The condition for the driving force 
to maintain the same direction as it has in GR is exactly 
eqn.( 151. It is interesting to see that in parameterization 
A Hs and C can have counteracting effects on the driv- 
ing term. Qualitatively, a negative ^-value with large 
magnitude enables one to weaken fi p considerably whilst 
maintaining growth during a matter-dominated epoch. If 
parameterization B is adopted this effect does not exist 
because the modification to the Poisson equation has no 
influence on the sign of the driving term. 

The friction term has a somewhat simpler behaviour, 
as it is unaffected by ^. In parameterization B the fric- 
tion is unchanged from GR, but in parameterization A 
fi = fip = fig has the power to enhance or suppress fric- 
tion effects. A value of /i > 1 acts to increase friction, but 
simultaneously strengthens the driving term that drives 



perturbations to collapse. One expects that some degree 
of cancellation between these two effects may be possible, 
even without the extra freedom provided by C. 



B. Super horizon scales 

On large scales we ought not to make the approxima- 
tion Am ~ Sm , as the magnitude of the velocity potential 
is not negligible. We will adopt a different strategy, us- 
ing the linearized Friedmann equation to solve for $ as 
a proxy for 5m- 

We first adopt a phenomenological approach (in the 
spirit of parameterization B), and assume that Newton's 
constant is identically modified in all perturbed field 
equations. We then have (neglecting a small term pro- 
portional to fc^) : 



-•H($ + -H*) 



6 



-PMfJ'pSM 



flpSM 



(16) 



Differentiating, and using that Sm ~ 3<i> on large scales 
(from eqn.(l6])): 



(17) 



$+^(l + 3/ip)$ + -H* = 
Using the slip relation (eqn.([5 
$ + $^[3 + 3^jp-2C]=0 



(18) 



The power-law solutions are $ cx , with = and 
— (2C — 2 — 3/ip). On these scales Sm follows the be- 
haviour of $ up to a constant offset, which can be set 
to zero by initial conditions. So the dominant mode - 
constant potential outside the horizon and Sm frozen - 
is the same as in GR. However, the decaying mode is 
affected by the MGFs. For C = 0, increasing fj,p will 
result in faster decay. This seems somewhat counter- 
intuitive, since one would usually associate an increase 
in gravitational strength with reduced decay of overden- 
sities. However, since we are working on superhorizon 
scales gauge issues may invalidate our physical notions 
of gravitational growth and decay. 

In parameterization A the linearized Friedmann equa- 
tion is 



-6H($ + -H*) = kq^pmSm 
where 

'fip - 1 



^0 



with Q 



ftp 
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(19) 



2C 



n 



The derivation of the above expression is described in 
appendix |B) 



Note that /ip does not feature explicitly in eqn.(19l 



Newton's constant is not modified directly - instead 
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Given the restriction that fig niust be positive, eqn.( 15 ) 
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FIG. 1: Joint constraints on the slip parameter and /ip 
at z — 0, for parameterization A (black lines) and B (filled 
green). In both cases 68% and 95% contours are shown. 



one considers all possible additional terms that could 
appear in the linearized Einstein equations, which can 
be determined up to a dimensionless function of back- 
ground quantities . In the case of second-order metric- 
only theory the only terms that can be added to the 
linearized Friedmann and 'Oi' equations are proportional 
to $, which can be absorbed into an effective Newton's 
constant. 

On very large scales limA;_j.o(fc^^o) = 0, so we can ne- 
glect the extra term on the RHS of eq.(19). Repeating 



the steps we took for the parameterization B-like case we 
obtain p+ ~ 0, = jip (2^ — 5). The values of in 
the two parameterizations converge as one tends to the 
GR limit, as of course they must. 

However, the decaying mode is not hugely interesting 
as it is unobservable (unless there are some very radi- 
cal modifications to GR involved). The important result 
here is that the potential remains constant on superhori- 
zon scales, as usual. 



tells us to expect that C < 2 on subhorizon scales in pa- 
rameterization A. For C = we must have /ip = /is > 1/2 
for density perturbations to grow during the matter era, 
which corresponds to the approximate location of the 
near-vertical black contour. This contour is not the re- 
sult of any artificially-imposed boundary, but delineates 
a very sharp fall-off in the likelihood distribution for /ip 
(see figure 6 of |23J)- In contrast, in parameterization 
B any /ip > permits growing modes for ^ = 0, giv- 
ing rise to the more gradual fall-off shown by the shaded 
countours. 

However, eqn.(15) also implies that in parameteriza- 



tion B C < 1 is necessary for growth of CDM density per- 
turbations on subhorizon scales, which is contradicted by 
Figjl] This is not too surprising - we expect the simple 
bounds implied by our EdS example to be blurred by the 
complexities of a realistic cosmological model. If models 
in the region of parameter space C > 1 experience suffi- 
cient growth during radiation- and A-dominated epochs, 
or on scales outside the validity of eqn.(14), then they 



will not be excluded by an MCMC analysis. 

The degeneracy between /ip and Q is visible in both 
contour sets, but it is more pronounced in parameter- 
ization B. This is because in parameterization A the 
quadratic term in eqn.(13) makes it more difficult to ac- 



commodate the effects of a large /ip with a small C, or 
vice- versa. 



D. Other Growth Observables 

Integrated Sachs-Wolfe Effect 

A cosmological model governed by GR will not experi- 
ence an Integrated Sachs- Wolfe (ISW) effect [SU] during 
a matter-dominated epoch. The kernel of interest for the 
ISW effect is $ -|- \1/ (in the conformal Newtonian gauge) , 
which in GR is equal to 2$. From the standard Poisson 
equation one has, on scales well below the horizon: 



$ (X o^PmSm 



(20) 



C. Connection to Constraints 

To what extent are our results for an idealized EdS 
model borne out in the real universe? Figjl] shows the 
joint constraints on the MGFs //p and C obtained us- 
ing the following data sets: the 7-year WMAP CMB 
data [45I, the SDSS DR7 matter power spectrum [46] . 
a prior Hq = 73.8 ± 2.4 [47 , the BBN constraint f^b/i^ = 
0.022 ± 0.002 48J, and the Union2 Supernova la data 
|49j . The differences between the two sets of contours 
were discussed in detail in [23]. Here our main interest 
is the extent to which the constraints reflect the analytic 
solutions of the previous two subsections. 



Energy-momentum conservation gives pM oc a~^, whilst 
in GR 5m oc a, leaving $ with zero time-dependence. 
The situation only changes for z < 0.5 when A begins to 
suppress the growth rate of 5m [51j . 

In modified gravity this behaviour is affected in two 
ways: a non-trivial slip relation will cause the ISW ker- 
nel to differ from 2$, and the scaling of 5m with a will 
be altered. Then generically we expect a non-zero ISW 
effect, even during a matter-dominated phase of the uni- 
verse [5H [53] . We will calculate the contribution to the 
CMB temperature power spectrum of this effect for the 
metric theories discussed in §III[ 

We begin by using the slip relation (eqn.([5])) to express 
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the ISW kernel purely in terms of 



$ + * = i$ 



3- 2C- 



1 



- $ 



- 1 



(21) 



Using the Poisson equation to connect $ and Sm, we have 
(discarding the decaying mode): 



where M{k) 



2fc2 



n+ is given by eqn.(12l and density perturba- 
tions are normalised by their present values i.e. 
SM{k,r]) = SM.oikjrj^^ . Generally one expects the 
growth rate of density perturbations to be scale- 
dependent in a modified gravity theory, whereas in GR 
all linear subhorizon modes grow at the same rate. How- 
ever, we are modelling a region of theory space close to 
GR and hence we will assume negligible variation of the 
growth rate over the range of k relevant to observations 
of the ISW plateau. 
So the ISW kernel is: 



$ + * = 



M{k) 



(71+ - 2)r7("^-3) X 



6- n+ - 2(+—{n+ - 2) 

Ms 



(22) 



Note that in the GR limit n+ = 2 the above expression 
vanishes as expected. 

Next we need to compute the power spectrum of this 
modified-gravity induced ISW effect. The expression to 
be evaluated is: 



Ci 



dkk^P{k) 



5Mfi{k) 



(23) 



The temperature perturbation observed today, Qi{k, r]o) 
consists of three parts: 



e(fc,?7o) 



monopole term -|- dipole term 



(24) 



drje $(fc,77) -f ^'(fc,?]) ji[k{r]o - f])] 



Let us focus on the dominant contribution to the ISW 
power spectrum which comes from the ($-1-^')^ term, and 
ignore the cross terms with the monopole and dipole. The 
cross-terms should only yield small corrections because 
the monopole and dipole terms are evaluated at the time 
of last scattering, and hence affect different Z-values from 
the subsequent ISW effect. We will take the visibility 
function e^'^ to be a step function at recombination. Of 
course, if we are considering a truly EdS universe then 
there is no recombination event or time of last scattering, 
but this detail is irrelevant - we are only interested in 
modelling the real universe well after recombination. In 
the real universe recombination occurs sufficiently early 
that for practical purposes we can take r^rec ~ 0. 

A standard derivation relates the power spectrum of 
density fluctuations today to the power spectrum of the 



/j = 1.0 f = -0.25 

1^=2.0 f=0.0 

IJ=2.0 f = -0.25 

- - /j=3.0 f =0.0 




FIG. 2: Low-Z power spectrum of the ISW effect induced by 
a theory of modified gravity constructed in parameterization 
A (see 311] for details). The y-axis scaling is arbitrary. The 
region I ^10 is expected to be subject to corrections. 



primordial potential [ST. The only modification to this 
that occurs in our theory is a factor of /ip^, which arises 
when we relate density fluctuations to $ via the Poisson 
equation. However, this is cancelled by a factor of /ip in 
M{k). The scales of interest to us are sufficiently large 
that we can set the transfer function T{k) ^ 1. Then, 
for a Harrison-Zel'dovich spectrum we find 

= ^(n+-2)2 (^6-n+-2C+^(n+-2)^'4 



dk 



(25) 



where Sh is the amplitude of primordial perturbations 
at horizon-crossing during a single-field slow-roll inflation 
scenario. The superscript 'sq' on C; reminds us that we 
have only evaluated the ISW-squared term and not the 
cross-terms. 

A plot of this power spectrum for several combinations 
of /i and C in parameterization A is shown in Figj2] The 
normalization of the y-axis is arbitrary because we have 
not attempted an accurate calculation of the C;s; we are 
more interested in how the shape and amplitude of the 
power spectrum is affected by different parameters. Note 
that our subhorizon solution for $ is likely to become 
invalid for the largest scales, hence the region / < 10 in 
Figj2]is not fully accurate. 

Provided that fip is not unusually small, can have 
a significant effect on the amplitude of the power spec- 
trum. For example, compare the curves with {/ip,C} 
equal to {1, —0.25} and {2, 0} - a small change in ^ domi- 
nates over a large change in fip. The models {2.0, —0.25} 
and {3.0, 0.0} have the same value of n+, and hence the 



7 



same spectral shape. It may seem a little surprising that 
a model with parameters {1.0,-0.25} predicts a larger 
ISW effect than one with {2.0, -0.25} for I > 6; this is 
because a larger value of 7i+ shifts the dominant contri- 



bution to the integral in cqn.(25 1 to later times, therefore 



shifting the corresponding power spectrum left towards 
larger scales. 

In general, however, it is difficult to cleanly disentan- 
gle the effects of and ( on the ISW power spectrum, 
because they appear in a degenerate combination inside 



Growth Function 

The rate of growth of structure as a function of redshift 
is often quantified via the growth function: 



din A 



M 



din a 



(26) 



In GR f{z) is independent of the wavenumber k, but 
in a modified gravity scenario this is not generally the 
case. However, if we assume we are not dealing with very 
radical departures from GR then this scale-dependence 
is likely to be small over a restricted range of k. On the 
scales of region 3 (see the beginning of this section) where 
Am ~ <5m) the growth function is simply f{z) = 
where is given by eqn.( 12 1. 



In a fully general case one should allow the effec- 
tive fluid to have a time-varying equation of state, non- 
adiabatic perturbations and significant anisotropic stress. 
This would render our simple EdS model invalid by 
modifying the evolution of the cosmological background. 
To proceed with our analytic treatment we will sacri- 
fice some generality by setting the equation of state of 
the effective fluid to zero. It then contributes to the 
zeroth-order Friedmann equation and supports density 
and velocity perturbations, but has negligible pressure 
and anisotropic stress. 

This approximation is not unreasonable; for example, 
in Hofava-Lifshitz gravity [551 [55] a non-local Hamilto- 
nian constraint gives rise to an integration constant which 
may be interpreted as a pressureless fluid. The additions 
to the Einstein equations in linear Einstein-Aether the- 
ory also behave as a pressureless fluid during an EdS 
phase, as can be seen using eqns.(15) and (16) of (60j . 
The effective fluid of Eddington-Born-Infeld gravity can 
experience a CDM-like phase too [BT]. However, making 
this restriction does exclude some important classes of 
theories, such as f{R) gravity, except for special choices 
of /(i?). 

The second type of new term that arises when extra 
d.o.f. are included is the following combination of metric 
potentials: 



r = ^ ($ + -H^' 



(27) 



IV. THEORIES WITH ADDITIONAL DEGREES 
OF FREEDOM 

Many gravity theories introduce degrees of freedom 
(d.o.f.) other than perturbations of the metric - addi- 
tional scalar or vector flelds, second metrics, or any com- 
bination thereof For model- independent constraints 
to be genuinely feasible we need to construct a param- 
eterization that is able to accommodate such d.o.f. In 
this paper we will restrict ourselves to scalar d.o.f. only, 
which flnd widespread motivation from particle physics, 
brancworld models and string theory. 

In SjiT] the parameterization- A-based approach was to 
add to each linearized Einstein equation all possible 
terms that could appear in the context of a metric-based 
second-order theory, which amounted to terms in $ and 
$ multiplied by some function of background quantities 
(see appendix [b|. When we extend this framework to in- 
clude extra d.o.f. two new types of terms appear. Firstly 
we must allow for perturbations of the d.o.f themselves 
- for example, in a Brans-Dicke theory [551 one expects 
the perturbations of the scalar fleld Scf), 5(j) and (50. Such 
new scalars are awkward to work with directly without 
knowing their underlying equations of motion. Therefore 
we will adopt an alternative approach based on a scheme 
by Hu [21] and treat this flrst type of additional term as 
perturbations of an effective fluid. A similar approach 
was adopted in [SS] |57j . 



In this paper we are working in the conformal Newto- 
nian gauge, but the parameterization we are using is 
constructed using gauge-invariant quantities - see |62| 
for details. In the gauge-invariant formalism $ and ^ 
are replaced by the Bardeen potentials —^h and ^a- 
If the background equations are fixed to be ACDM, as 
in pil[ then the gauge-invariant version of F cannot ap- 
pear in the linearized field equations because it intro- 
duces higher-order time derivatives of the scale factor. 
However, when new d.o.f. are added these unwanted time 
derivatives can be cancelled by perturbation of the new 
scalar. This somewhat technical point was demonstrated 
explicitly for scalar-tensor theories in appendix B of [22j . 

Concretely, then, we write the Einstein equation in the 
form: 

G^, = STrGoa^T^. + a^f/^, (28) 

where the tensor U^i, contains all the non-standard terms 
arising from a theory of modified gravity. At the per- 
turbed level the components of 5U^u are then given by 
(in Fourier space): 

L/a = A^k'^f^ + Fok^T + na'^pEdE 

Ue = Bok^ + lokT + ko^peOe 

Up = Cofc^$ + Cifc<i> + Jofc^F + Jifcf 

Us = Do^ + + KoT + (29) 
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where 

Ua = -a^SU^, 
Up = a^SW, 



D,,U^ = a^iSU; - hu^S)) (30) 



The coefficients Aq, . . . Ki above are functions of back- 
ground quantities, dependencies which we wih suppress 
to avoid cluttered expressions. A subscript E denotes 
the quantities relating to the effective fluid. The expres- 
sions above contain fewer terms than the corresponding 
ones in 1221 : the additional metric terms present in that 
paper are those needed to form a gauge-invariant com- 
bination with perturbations of the new scalar. Here we 
have kept those terms folded into the effective fluid so 
that the correspondence with SUf^i, in a general gauge is 
clearer. 

The system to be solved then comprises of the six vari- 
ables {^,T,Se,Oe,Sm, Sm} and six dynamical equations: 
the two spatial components of the linearized Einstein 
equation and two fluid conservation equations for each 
of CDM and the effective fluid (equivalent to Bianchi 
identities). Note that because SU^^ contains additional 
metric perturbations the conservation equations of the ef- 
fective fluid will contain some non-standard terms. The 
full system of equations is displayed in appendix |Aj 

r can be eliminated from the two spatial Einstein equa- 
tions to give a second-order equation in terms of $: 



a 
a 



a 



Wi 
' W~i 

W2- 



kZi 
a 



(31) 



Wi 



W2Z1 , WiZ: 



k- 



where Hk = H/k and 

a = (^Di- ( Ji - 6) -I- i^i (dHk ~Ci + ^^ 

Wi = Ki (^Uk {Ji - 6) (^Ko + 

W2 = {Ji - 6){Do - I) + Ki{2 - Co) 

Z2 = {2- Co) (:^^ - + {Do - 1) ( mk - Ci 



2 



We have assumed that time derivatives of the MGFs are 
negligible, consistent with our treatment of metric-only 
theories. If we set the MGFs to zero the third term van- 
ishes and we recover the GR result $ = constant. Modi- 
fications to the Poisson equation mean that the solution 
for $ is not as easily translated into a solution for 5m as 
it is in GR. One route is to substitute the solution for $ 



both of these solutions in eqn. 
an explicit calculation. 



Q - see appendix [C] for 



We wish to consider eqn.(|31| on superhorizon and sub- 
horizon scales, 



as we did in §III| However, there is some 
uncertainty involved in taking this limit without know- 
ing specifically the functional forms hiding behind the 
MGFs. We will assume that any scale-dependence ap- 
pears relative to some preferred scale of a given theory, 
i.e. as a function of (fc/fc*). 



A. Subhorizon Scales 

Under the assumptions stated above, and retaining 
only the dominant terms when expanded in powers of 
(k/H), on subhorizon scales eqn. (31) reduces to: 



{2Di - Ci 



{6-J1+2K1) 



-k<i> + 



{6-J1 + 2K1] 



fc^$ = 



(32) 



With the change of variable x = krj the above equation 
can be rewritten (using primes to denote derivatives with 
respect to x): 



7 

where P = 



7 

{2Di -Ci + 



$ 

2i^o-Jo), 7 



6- Ji 



(33) 
2Ki. 



Without taking specific forms of the MGFs we cannot 
say whether the coefficients of the second and third terms 
above are positive or negative, only that sufficiently large 
modifications to gravity have the power to flip their signs 
(although the magnitude of such non-GR terms is ex- 
pected to be small in the domain of validity of this equa- 
tion). However, one might expect (2Z?o — Co)/7 > and 
/3/7 > so that changes in $ damp out (and therefore 
return to the GR-like situation) rather than grow. 

Eqn.([33| has the form of a simple mechanical oscilla- 
tor, and will display the usual phenomena of ringing or 
over-damping depending on the values of the coefficients. 
Specifically, its behaviour will depend on the value of 
/3^/7 — SDq + 4Co, with negative values of this quantity 
leading to damped oscillations in $(a;) and positive val- 
ues leading to exponentially growing and decaying solu- 
tions. Let us apply the single boundary condition that 
the potential is constant on superhorizon scales; we will 
see shortly that this is likely to remain true in theories 
with extra d.o.f. We can then determine the subhorizon 
solution up to an overall constant: 



$(a;) cx e' 



where m± 




(34) 



into either of eqns.(Al) or (A2) and solve for F, then use 



At first sight this exponential solution might seem to be 
a cause for concern, as we would normally expect the 
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growth of <I> or 5m on small scales to follow power-law 
behaviour. The unusual solution above has arisen be- 



cause the modifications in eqns.(29) introduce factors of 



k that dominate over the usual GR terms on small scales. 
In the toy model considered here we took the MGFs to 
be dimensionless functions of order one; without such as- 
sumptions it is difficult to make any general statements 
about the growth of perturbations, because we do not 
know how to assess the relative importance of two fac- 
tors such as Jo and 3k/T-l. Of course, if we know the 
specific functional forms hidden behind the MGFs then 
such assumptions are not necessary, but then we would 
not be pursuing a model-independent approach. 

A more realistic situation would be to take the MGFs 
to be much smaller in magnitude than the GR terms. 
However, our intention here is to assess qualitatively the 
effects that parameterized systems of modifications to 
gravity have on the growth of structure - a task which 
becomes difficult when the parameters are taken to be 
vanishingly small. We will therefore maintain the sim- 
ple assumptions described above, remembering that in 
a more realistic model the effects described here would 
only be manifest as small distortions of a predominantly 
GR-controUed universe. 

Alternatively, we could turn this problem on its head. 
One reason theories with additional d.o.f. are difficult 
to work with is because we lose the ability to derive a 
hierarchy of constraint equations between the MGFs, as 
we did in the purely metric case (these arise from the 
Bianchi identities - see appendix [b|. Can we use the 
growth of $(a;) to infer some replacement constraints? 

To clarify, we wish to restore power-law behaviour 
for ^(x). By expanding eqn.(31) in powers of k/T-L we 
can find the conditions necessary to remove the domi- 
nant terms that are causing the exponential solution of 
eqn.(35). We find these to be: 



2Do - Co = 0, 2Di - Ci = 0, 2Ko - Jo = (35) 
If the above conditions are satisfied on small scales then 



eqn.(31 1 reduces to: 



I -Do 



12 



$(l-i?o) 



6 



- 1 



= 
(36) 



results: 
dMia) = Pk 



(Do - 1) 



{q+ + 2)iq++3) 
Pk"" 16(1- Do){l-flE) 



(38) 



{q+ + 2){q+ + 3) 



1 

MP 



a^+i (39) 



In the GR limit Qe ^ 0, Dq ^ (which sends q+ 0) 
we recover that Sm scales with a. 



It is interesting to see that the conditions in eqn. ( 35 ) 
are satisfied in parameterization A of Sjll] in the limit 
fc — > oo (the last condition trivially so) . The same is true 
in parameterization B, where Di — 0. If one wished to 
implement this gravitational framework into a realistic 
cosmological model, without necessarily taking the ab- 
solute magnitude of the MGFs to be very small, then 
these are the constraints that must be satisfied to give 
a reasonable degree of structure formation. Of course, 
this does not prevent the model from being ruled out 
by other observables such as the ISW effect. Having re- 
stored power-law growth, the effects of this theory on 
the ISW effect and growth function are expected to be 
qualitatively similar to those in §111 D| 



B. Superhorizon scales 



On superhorizon scales eqn. (31) simplifies to 



18Ki 



Ko- 



(40) 



iJi-l2){Do-l) + Ki (2 -Co 



= 



This can be reduced to an oscillator equation in krj, in 
complete analogy to the subhorizon case. However, in 
the superhorizon limit A: — > oo these oscillations become 
infinitely slow in fcry; effectively, we have that <I> is a con- 
stant. This matches the GR and metric-only cases. 



CONSERVED SUPERHORIZON 
PERTURBATIONS? 



The solutions are then power laws in a (or rj), as desired: 
$(a) = Pa'^ + Qa^ 

In appendix [C] we convert this solution for $(a) into a 
solution for Sm and Se- Here we shall simply state the 



Any relativistic theory of gravity in which energy- 
momentum is covariantly conserved allows the definition 
of a perturbation, Z, that is conserved on superhorizon 
scales in the absence of non-adiabatic perturbations 
I65j . In the literature Z is more commonly denoted as C, 
but the choice of notation for one of the MGFs in this 
paper and previous work prevents us from reusing that 
letter. Z is identified with the curvature perturbation 
on uniform-expansion hypersurfaces in a homogeneous 
and isotropic spacetime, which in GR coincide with hy- 
persurfaces of constant energy density. For zeroth-order 
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Einstein equations of the form: 



- " f 



-9o 



the conserved perturbation is [64] : 
fo 



(41) 



(42) 



It would be interesting to know how Z behaves in the 
metric theories considered in 31111 Since the new non-GR 



terms do not really originate from perturbations to a 
fluid, it is not immediately obvious whether they will 
be equivalent to adiabatic or non-adiabatic pressure 
perturbations. Therefore the conservation of Z does not 
necessarily follow. 



The perturbations SPm and SY have been decomposed in 
a manner analogous to eq.(44). Fp^^x represents a pos- 



sible entropy perturbation between CDM and the modi- 
fied sector, which can be non-zero even if each component 
does not support entropy perturbations by itself. Pertur- 
bations of the background quantities X and Y correspond 
to components of the tensor [7^,^: a^SX = —o^SUq and 
a^SY — a?6Ull?>. Then the effective non-adiabatic pres- 
sure perturbation of the modified sector is: 



5Y 



5Ul 
3 



(48) 



Our toy model of an EdS universe contains only cold dark 
matter, so 5Pm — SPnad — c^^*^^ — 0. 



Using the linearly perturbed versions of eqns.(41 1, one 



can derive an equation for the evolution of the metric 
potentials [64] : 



mn -U-W 



2n^ - nn 



a 
~2 



(43) 



where the perturbation 5g has been decomposed into 
parts equivalent to adiabatic and non-adiabatic pressure 
perturbations: 



Jo 



inad 



(44) 



The time derivative of Z is related to the quantity on the 
LHS of eqn.(43|. Hence the rate of change of Z is found 
to be: 



Sgnad 



(45) 



Uii ~ Y we can read off: 

/o = + X 

5o = k(pm + Pm 



Comparing eqns.(28) and (41) and defining Uoo — X, 



X + Y 



(46) 



Rearranging eq.(44) and substituting in our expressions 
for / a 

Sgnad — 



for / and g, one finds 
kpmX 



KPM +X 



A. Metric-Only Theories 

For a purely metric theory with a ACDM-like back- 
ground X and Y are zero or equivalent to a cosmological 
constant (i.e. X + Y = 0), so there can be no entropy 
perturbations between CDM and the modifications (since 
Cs —0). The only possible source for non-conservation 
of Z would be from non-adiabatic perturbations within 
the modified sector itself, SY^ad- Eqn.(45) becomes: 



7iad 



(47) 



where 
TpMX = 

^2(Af) _ 



SpM 

Pm 

Pm 

Pm 

Y 

~X 



5X_ 
X 



n 



2(^-^2) 3 

n 



(49) 



Non-adiabatic perturbations within the modified sector 
would amount to fluctuations about uje = —1, which 
could lead to a situation equivalent to a phantom field. 
Whilst a phantom equation of state is permitted by cur- 
rent data [5SH55] , direct phantom scalar field models are 
plagued by severe difficulties because they lead to an un- 
stable vacuum state [69J and favour an anisotropic uni- 
verse }70j . However, it is known that theories such as 
scalar-tensor gravity, /(i?) gravity and some Lorentz- 
violating models can cause phases of an effective uje < — 1 
without introducing a phantom field per se |71H75] . 

Fortunately the question of whether non-adiabatic per- 
turbations lead to < —1 turns out to be a moot point 
for the theories considered in ^[TT] and pilj This is be- 
cause the effective pressure perturbation in eqn.(48) van- 
ishes on very large scales. The functions Co and Ci are 
related to pp and C by a set of constraint equations, dis- 
played in appendix [B] There we also show that under the 
assumptions mades in this paper limfc_j.o (fc^Co) = and 
limfc_>.o (fcCi) = in both of the parameterizations de- 
scribed in [jlT] . These results agree with the conclusions 

of 
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B. Theories with Additional Degrees of Freedom 

For theories with extra scalar degrees of freedom the 
situation is different. Firstly, the final term in SYnad 
does not vanish, as the modifications to the background 

'2( J^) 

equations lead to 7^ 0. Secondly, we cannot derive a 
system of constraint equations on the functions Aq . . . Ki 
like those in Table |T] because the new scalar now acts 
as a source in the Bianchi identities (see appendix |b] for 
details). Thirdly, we can now have entropy perturbations 



between matter and the modified sector, as the coefficient 



of ^pmX in eqn.(47) no longer vanishes. 



Maintaining the approach of pV[ we keep the metric 
components oiSUf^i, distinct, but treat the perturbations 
of the new d.o.f. as an effective fluid. We previously 
set the equation of state of this effective fluid to zero 
in order to preserve the EdS nature of our toy model; 
let us restore the general case for the present. The full 
expression for Z then becomes: 



n 



X 



X 



\ UPM + X 



X 



na^PE i^E 



X 



1 . 



npM 



X 



(50) 



If the MGFs do not contain inverse powers of k then 
the second line vanishes on very large scales, leaving only 
perturbations to the effective fluid (although a similar as- 
sumption was used in \1Y to make an analytic solution 
achievable, it does not necessarily have to hold true for 
all theories). The first set of round brackets represents 
entropy perturbations between CDM and the modified 
sector, and the second term is akin to non-adiabatic per- 
turbations within the modified sector itself. However, 



\ IV in which the effective fluid was assumed to be pres- 



note that 



is not necessarily equal to the sound 



speed of the effective fluid, since X and Y may con- 
tain terms constructed from the metric. This feature 
distinguishes a dark energy model like quintessence from 
a modified gravity model such as a scalar-tensor gravity. 
In the latter the scalar field is non-trivially coupled to 
the metric and hence X and Y contain more terms than 
just the energy density and pressure of a scalar field. A 
quintessence-like case can be recovered from the above 
expression by setting the MGFs to zero and X = KpE, 
Y ~ K Pe , which gives: 



Z = - 



where 



m.e 



n 



/9_e(1 +UJE)dt''^^TM,E 



1 + {1+uje) 



E ,nad 



3M 



l+UE 



(51) 



(52) 



If the quintessence field supports only adiabatic pertur- 
bations {5PE,n&d = 0) then Tm,e can be set to zero 
through choice of adiabatic initial conditions, leaving Z 
conserved. 

As a second example, consider the theory treated in 



sureless. Eqn.(50) reduces to: 



X 



(h-H^^ 2 (^KpM + X 



cl'-^^ [pmSm + PeSe] 



(53) 



The square brackets can be set to zero at through a choice 
of isocurvature initial conditions, so Z remains conserved 
a uje is constant. However, as one wishes the effects of 
modified gravity to become apparent at late times in the 
universe an evolving ue would be more desirable, for 
which Z would not be conserved. 



VI. CONCLUSIONS 

Awareness of a forthcoming 'data deluge' from cur- 
rent and future cosmological experiments has led to an 
interest in model-independent approaches to constrain- 
ing modified gravity. By exploiting the generic features 
of current models, these parameterized systems seek to 
constrain large regions of theory space simultaneously. 
Whilst there has been much effort made to constrain 
these parameterizations with the data [25l EH ES] , there 
has been relatively little investigation into the corre- 
sponding theoretical description: how does the param- 
eterized system of perturbation equations evolve? (Al- 
though see [21]; general scalar field-type models are 
treated in [H [771 [75] ) • 

We have attempted to answer this question within the 
simplified setting of an Einstein-de Sitter universe, for 
two classes of theories: 1) those for which the degrees 
of freedom are the metric and matter perturbations, and 
2) theories which explicitly introduce additional scalar 
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degrees of freedom. In both cases we have found that 
the curvature perturbation remains constant on the very 
largest scales (well above our observable horizon). 

On subhorizon scales we found that in the metric-only 
case perturbations grow in familiar power-law fashion, 
but the exponents are modified from their General Rel- 
ativistic values. The two modified gravity functions fip 
(controlling the effective Newton's constant) and ( (con- 
trolling the dominant component of the gravitational 
slip) are partially degenerate in their effects, making 
them difficult to disentangle. The impact of C becomes 
subdominant if the effective Newton's constant is weak- 
ened. Modifications to the evolution of $ lead to an 
induced ISW effect and modification to the growth rate. 

Theories with additional scalar degrees of freedom are 
considerably harder to study, as the number of undeter- 
mined functions is much larger in this case. In §IV| we 
considered phases during which the new terms in the Ein- 
stein equations behave as an effective fluid with a negligi- 
ble equation of state. We found that the modified terms 
dominate the evolution equations, leading to damped os- 
cillatory or exponential behaviour. From a study of these 
equations we have found three relations (eqns.( [35| )) be- 
tween the modified terms that, if satisfied, will restore 
power-law behaviour. This analysis docs not apply to all 
modified gravity theories (when w^; is not small), but it is 
relevant to (for example) EdS regimes of linear Einstein- 
Aether theory, EBI gravity and Hofava-Lifshitz gravity. 

Constraint equations such as this are desirable be- 
cause the number of free functions required to param- 
eterize common gravitational theories increases rapidly 
when new degrees of freedom are included. This prolifer- 
ation could reduce our ability to constrain such param- 
eterized frameworks satisfactorily. One can reduce this 
freedom somewhat by treating the new scalar degrees of 
freedom as an effective fluid; indeed this (or an equivalent 
approach) is necessary if the parameterization is to cap- 
ture theories which introduce more than two new scalar 
degrees of freedom. The metric perturbations are also a 
considerable source of freedom, as displayed in eqns.(p9 1. 



This splitting of the modifications into metric parts and 
effective fluid parts is the key to distinguishing between 
closely related models of dark energy and modified grav- 
ity, such as quintessence and scalar-tensor theories. 

Under the assumptions made in this paper the met- 
ric terms become irrelevant on ultra-large scales. In a 
metric-only theory this leaves the superhorizon pertur- 
bations Z (more commonly denoted by by C.) conserved, 
but in class 2) theories the possible non-conservation of 
Z depends on the equation of state of the effective fluid. 

The EdS model considered in this paper is useful in 
obtaining a qualitative understanding of how parame- 
terized gravity might effect the matter-dominated phase 
of our universe. To obtain more quantitive predictions 
this must be embedded in a more complex cosmological 
model, which is likely to be achievable only numerically. 
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Appendix A: System of Equations for Theories with 
Extra D.o.F. 



This appendix displays the system of six equations that 



is solved in ^IV These are: the two spatial components 
of the Einstein equation (longitudinal and transverse tr- 
acless), the fluid conservation equations for cold dark 
matter (the two components of 6 {W ^Tj^) = 0) and the 
two Bianchi identities S {W/^iUj^) = 0. The non-standard 
terms containing MGFs arise from the parameterization 
laid out in eqns.(28) and (29), which treats the additional 



scalar degrees of freedom as a pressureless fluid (denoted 
by a subscript E). This reduces to parameterization A 
used in §III| for purely metric theories. 

The variable ^E* in these equations can be eliminated in 
favour of F using the definition F = ($ -I- 'H'^)/k. 



6fcF -I- 12-HfcF -I- 2k^ ($-*)+ 6(-H - H^)* = Cofc^$ 4- Cik^ + J^k^T + J^kV (Al) 

= L>o* + + A'or-I- (A2) 
k k 

Sm = -k^eM + 3^ (A3) 

eM = -neM + ^ (A4) 

ku^PeSe = Ka'^pE (-k^OE + 3$) - fc^^ {UAq -f kB^ + HC'o) + fc$ {kAo + UCi) + k^T {UFq + klo + HJq) 

+ktikFo+HJi) (A5) 

Ka^PE0E = «a'pB i-neE + *) - - Ico + Ido] - fc$ (b„ - + 



-k'F I 2^/0 - :^ Jo + Iko] - kt (lo - + ] (A6) 
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Constraint equation 


1 


io + HAo + kBo + HCo = 


2 


Ao+mCi =0 


3 


Bo + 2HBo - \kCo + ^kDo = 


4 


Bo - sCi + ffi = 



TABLE L Constraint equations for a metric theory in param- 
eterization A, specified by eqnsJBll. 



Appendix B: Conservation of Z in Metric-Only 
Theories 



In this appendix we demonstrate that any effective 
non-adiabatic perturbations that might prevent Z from 
being conserved vanish on large scales for the metric- 
only theories considered in and . First we must 
introduce some more detail about the parameterization 
underlying eqns.Q and ([5|. 



Consider a modified Einstein equation of the form ( 28 1 . 
In |22] we demonstrated that in parameterization A (in 
which fig = fip) a theory containing up to second-order 
time derivatives corresponds to the following forms for 



the components of SU^i, (see eqns.(30l for definitions of 
quantities on the left-hand side): 



Ua 
Up 



Bok^ 

Cok^^ + Cik^ 
Di ■ 



k 



(Bl) 



where the coefficients Ao,...Di are functions of back- 
ground quantities, dependencies which we will suppress 
to avoid cluttered expressions. $ is a gauge-invariant 
perturbation variable, equivalent to the Bardeen variable 
— ^'ij, which reduces to $ in the conformal Newtonian 
gauge. 

The two components of S (V^?7/f ) — result in equa- 
tions containing $ and its derivatives. To avoid contra- 
dicting the solution for $ dictated by the Einstein equa- 
tions these expressions must vanish identically. Setting 
the coefficients of each metric perturbation to zero leads 
to the four constraint equations listed in Table [ij This 
enables us to reduce the six free functions in eqns.( [BT| ) 
to just two. Note that terms in 4' (equivalent to the 
Bardeen variable $^ ) are not permitted to appear in SUq 
and SU^, because they contain second-order time deriva- 
tives, which would lead to third-order Bianchi identities. 
Constraint equations similar to those in Table |l] prevents 
\1/ from featuring in SU^ and SU^ if it is not present in 
6U^ and SU°. 

We will choose one of our free functions to be Dq - 
this corresponds to ^ in eqn.(|5|. The second free function 
will be the combination that appears when we form the 





Constraint equation 


1 


Ao + nAo + kBo + nco = 


2 


Ai + HAi + kAo + kBi + HCi = 


3 


kAi + HC2 = 


4 


Bo + 2HBo - ffcCo + pDo = 


5 


Bl + 2UBi + kBo - |fcCi = 


6 


Bl - IC2 = 



TABLE II: Table of the constraint equations for a metric the- 
ory in parameterization B, specified by eqns.(|B4|. 



Poisson equation: 
1 



1 + i (Ao + S-HkBo) 



(B2) 



where Hk — H/k. In terms of these two MGFs the other 
coefficient functions arc: 



Ao 



Bn = 



Co = 



MP - 1 



1 



MP / V Q 
2nk ~ c) 

3Q 

2 /(mp-I) 



2C 



fj,p 



(Mi 



MP 



Ci = 



Di = 




2C 



1 



2 

9WI 



2C 



J_ / mp- 1' 
T-Lk \ MP 
\_ /mp-1' 
T-LkX MP 

g = + ^ 

o 

The second expression for Co applies only in an EdS 



where 



-n. 



(B3) 



universe. 



From eqn.(49l we see that the quantity of 
interest for evaluating Z is 5Ul. In fact we only need to 
know Co, since we showed in §111 B that the potential is 
constant on large scales. Under the assumption tha t ^p 
and C do not contain inverse powers of fc, eqns.(Bl) and 
(B3) imply that Mmk^o 5Ul — 0. Hence Z is conserved 



on superhorizon scales in parameterization A. 

We can repeat this calculation for parameterization B, 
in which /i^. = 1 but /ip remains a free function. This 
corresponds to a SU^j^i, tensor of the following form (see 
[22]): 



Ue 
Up 



Aok"^^ + Aik^ 

Bok^ + Bi^ 

Cofc2$ + Cifcl> + Ca* 



(B4) 
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The constraint equations for this theory are hsted in Ta- 
ble [n] In terms of the MGFs fi p and ( the coefhcients 
are: 



Al = 



Bo = 



MP / V 
-GHk / /ip - 1 



V MP 



i 

3Q 

2 f^p-V 



{i + ml) \ MP y 



2 /(mp-1) 



Q V MP 



C)|^H^-H§ U2C 



^ 2m / (/ip~i) _ 



V MP 

.4mf(-^)I^^ 

2-Hfc /(mp - 1 



h4-Hfc 



MP 



-c 



Co 



MP / [Q-lH^y 
6 / /^p — l' 



(1 + 3^2) V MP 



(B5) 



As we found for the parameterization A case, 
Umfe_i.o 5Ul — 0. Hence Z is conserved in parameteri- 
zation B also. 



Appendix C: 5m and 5^ in Theories with Extra 
D.o.F. 



Eliminating V from these: 



18 



7 ( ^ 



Al 



7 fx 



X Ki \2 



(CI) 

= 



where 7 = 6 — Ji -f 2A'i. Keeping only the growing mode, 
we write the solution for $ as: 

$(2:) = i?(/fc)x« (C2) 
We have absorbed the fc-dependence into the prefactor. 



R{k) = P/k'' where P is a constant (see eqn.(37)), and 



dropped the superscript -I- on q to avoid cluttered expres- 
sions. 



Taking the subhorizon limit x — > 00, eqn.(Cl) gives us 
the following solution for F: 



r{x) = 2R{k) ( 



l-Do]x' 



(C3) 



Next we recast the fluid conservation equations for CDM 
in dimensionless format and combine them in a manner 
analogous to eqns.(|6])-(|9|, yielding: 



3$" 



/ 6 X 



2 



This has the solution: 



bu^x) = ci + -+R{k)x 



{Do - 
iq + 2){q + 3) 



(C4) 



(C5) 



where ci and C2 are integration constants. This solu- 
tion has the desired behaviour that when we set Do — 



(which sets 9 = 0, see eqn.(37l) we recover that 6m grows 
as Sm oc a;^ cx a, as occurs in GR-controUed matter- 
dominated epoch. Finally we use the (dimensionless) 
Poisson equation to relate the solutions for $(a;) and 
(5a/ (a;) to (5b (x): 



We wish to connect the subhorizon solution for $ in 



eqn.(37) to the density perturbations of CDM and the 



effective fluid. First we cast eqns.(Al) and (A2| in di 



mensionless format by using the substitution x — krj (and 
correspondingly V. ~ 2/x, H = — 2/a;^). We also apply 
the conditions necessary for power-law growth, given in 
eqns.(35). The spatial components of the Einstein equa- 



tions become (where primes denote derivatives with re- 
spect to x): 



6 
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X - 2Di 



r'{6 - Ji) + r - X - 2Ko 

+ 2$(i-i:»o) = 

r'Ki + r (xo + 1) + {di - I) - $(1 ~Do)^o 



nE5E) + r^lp[Fo + -Io] (C6) 



-2$ = 1^ (1]m<5m 



In our toy EdS universe, containing only CDM and the 
effective fluid, ^Im and Qe are constants. Retaining just 
the dominant growing modes for both Sm{x) and 5e{x), 
the above equation gives us the result: 



dE(x) 



R{k) 
6Qe 



6(1 - Do){l - Qe) }_ 
((7 + 2)(g + 3) iip 



(C7) 



As expected, (5m (a;) and bE{x) grow at the same rate. 
Note that, as in GR, their evolution differs from that of 
$(a;) by a single factor of a. 



15 



A. Nicolis, R. Rattazzi, and E. Trincherini, Phys. Rev. D 
79, 064036 (2009), 0811.2197. 

V. A. Rubakov, Physics Uspekhi 44, 871 (2001), 
arXiv:hep-ph/0104152.' 

R. Maartens and K. Koyama, Living Reviews in Relativ- 
ity 13, 5 (2010), 1004.3962. 

T. Clifton, P. G. Ferreira, A. Padilla, and C. Skordis, 
ArXiv e-prints (2011), 1106.2476. 

S. Nojiri and S. D. Odintsov, Physics Reports 505, 59 
(2011), 1011.0544. 

C. M. Will, Astrophys. J. 163, 611 (1971). 

K. S. Thorne and C. M. Will, Astrophys. J. 163, 595 

(1971) . 

C. M. Will and K. Nordtvedt, Jr., Astrophys. J. 177, 
757 (1972). 

C. M. Will, ArXiv e-prints (2011), 1102.5192. 
L. Landau, Soviet J. Phys (JETP) 116 (1946). 
P. J. E. Peebles, ApJL 263, LI (1982). 

P. J. E. Peebles, The Large-Scale Structure of the Uni- 
verse (Princeton University Press, 1980). 
G. R. Blumenthal, S. M. Faber, J. R. Primack, and M. J. 
Rees, Nature (London) 311, 517 (1984). 
P. Meszaros, A&A 37, 225 (1974). 

S. Nesseris, C. Blake, T. Davis, and D. Parkinson, JCAP 
7, 37 (2011), 1107.3659. 

G. Bambi, M. Giannotti, and F. L. Villante, Phys. Rev. 

D 71, 123524 (2005), arXiv:astro-ph/0503502. 

E. Galabrese, D. Huterer, E. V. Linder, A. Melchiorri, 

and L. Pagano, Phys. Rev. D 83, 123504 (2011), 

1103.4132. 

G. L. Reichardt, R. de Putter, O. Zahn, and Z. Hou, 
ArXiv e-prints (2011), 1110.5328. 

P. Brax and A.-C. Davis, ArXiv e-prints (2011), 
1109.5862. 

G. Skordis, Phys. Rev. D 79, 123527 (2009), 0806.1238. 
P. G. Ferreira and G. Skordis, Phys. Rev. D 81, 104020 
(2010), 1003.4231. 

T. Baker, P. G. Ferreira, G. Skordis, and J. Zuntz, ArXiv 
e-prints (2011), 1107.0491. 

J. Zuntz, T. Baker, P. Ferreira, and G. Skordis, ArXiv 
e-prints (2011), 1110.3830. 

W. Hu, Astrophys. J. 506, 485 (1998), astro-ph/9801234. 
R. Bean and M. Tangmatitham, Phys. Rev. D81, 083534 
(2010), 1002.4197. 

S. F. Daniel and E. V. Linder, Phys. Rev. D 82, 103523 
(2010). 

L. Pogosian, A. Silvestri, K. Koyama, and G.-B. Zhao, 

Phys. Rev. D 81, 104023 (2010), 1002.2382. 

A. Hojjati, G.-B. Zhao, L. Pogosian, A. Silvestri, 

R. Grittenden, and K. Koyama, ArXiv e-prints (2011), 

1111.3960. 

E. Bertschinger and P. Zukin, Phys. Rev. D 78, 024015 
(2008), 0801.2431. 

W. Hu and I. Sawicki, Phys. Rev. D76, 104043 (2007), 
0708.1190. 

L. Pogosian, K. Koyama, P. Ferreira, and G. Zhao (2011), 
private communication. 

D. Lovelock, Journal of Mathematical Physics 13, 874 

(1972) , URL http: //link . aip . org /link/?JMP/13/87471 

D. Lovelock, Journal of Mathematical Physics 12, 498 



(1971), URL |http : //link . aip . org/link/? JMP/12/498/ 

I — n 

S. Deser and R. P. Woodard, Physical Review Letters 
99, 111301 (2007), 0706.2151. 

C. Deffayet and R. P. Woodard, JGAP 8, 23 (2009), 
0904.0961. 

Bruno and Zumino, Physics Reports 137, 109 (1986), 
JSSN 0370-1573, URL http : //www . scien cedirect . com/| 
"science/article/pii/037015738690076l| 

L. Amendola, M. Kunz, and D. Sapone, JGAP 4, 13 



[34 
[35 
[36 



L. 

[37 
[38; 
[39 

[4o; 

[41 

[42 

[43; 

[44 

[45' 
[46; 

[47 

[48 

[49; 

[50 

[51 



[52 

[53; 



(2008), 0704.2421. 

G. Schimd, J.-P. Uzan, and A. Riazuelo, Phys. Rev. D 
71, 083512 (2005), arXiv:astro-ph/0412120. 
L. Pogosian and A. Silvestri, Phys. Rev. D 77, 023503 
(2008), 0709.0296. 

K. Koyama and R. Maartens, JGAP 1, 16 (2006), 
arXiv:astro-ph/0511634. 

G.-B. Zhao, L. Pogosian, A. Silvestri, and J. Zylberberg, 
Phys. Rev. D 79, 083513 (2009), 0809.3791. 

A. Hojjati, L. Pogosian, and G.-B. Zhao, ArXiv e-prints 
(2011), 1106.4543. 

J. Dossett, M. Ishak, and J. Moldenhauer, ArXiv e-prints 
(2011), 1109.4583. 

T. P. Sotiriou and V. Faraoni, Reviews of Modern Physics 
82, 451 (2010), 0805.1726. 

D. Larson et al., ApJS 192, 16 (2011), arXiv:1001.4635. 

B. A. Reid et al., MNRAS 404, 60 (2010), 
arXiv:0907.1659. 

A. G. Riess et al., Astrophys. J. 730, 119 (2011), 
arXiv: 1103.2976. 

J. Hamann, J. Lesgourgues, and G. Mangano, JGAP 3, 
4 (2008), arXiv:0712.2826. 

R. AmanuUah et al., Astrophys. J. 716, 712 (2010), 
arXiv: 1004. 1711. 

R. K. Sachs and A. M. Wolfe, Astrophys. J. 147, 73 
(1967). 

R. G. Grittenden and N. Turok, Phys. Rev. Lett. 76, 



575 (1996), URL http : //link. aps . org/doi/10 . 1103/ 



PhysRevLett ■76.575( 

TT Zhang, Phys. Rev. D 73, 123504 (2006), URL |http : | 



//link . aps . org/doi/10 . 1 103/PhysRevD . 73 . 123504 

T. Giannantonio, M. Martinclli, A. Silvestri, and A. Mel- 
chiorri, JGAP 4, 30 (2010), 0909.2045. 
[54] S. Dodelson, Modern Cosmology (Academic Press, 2003). 
[55] C. Brans and R. H. Dicke, Phys. Rev. 124, 925 (1961). 
[56] T. Koivisto and D. F. Mota, Phys. Rev. D 73, 083502 

(2006), arXiv;astro-ph/0512135. 
[57] D. F. Mota, J. R. Kristiansen, T. Koivisto, and N. E. 

Groeneboom, MNRAS 382, 793 (2007), 0708.0830. 
[58] P. Horava, Phys. Rev. Lett. 102, 161301 (2009). 
[59] P. Horava, Phys. Rev. D 79, 084008 (2009). 
[60] T. G. Zlosnik, P. G. Ferreira, and G. D. Starkman, Phys. 

Rev. D 77, 084010 (2008), 0711.0520. 
[61] M. Baiiados, P. G. Ferreira, and G. Skordis, Phys. Rev. 

D 79, 063511 (2009). 
[62] G. Skordis, Phys. Rev. D 79, 123527 (2009). 
[63] D. Wands, K. A. Malik, D. H. Lyth, and A. R. Lid- 
die, Phys. Rev. D 62, 043527 (2000), arXiv:astro- 
ph/0003278. 

[64] A. Gardoso and D. Wands, Phys. Rev. D 77, 123538 
(2008), 0801.1667. 



16 



[65] E. Bertschinger, Astrophys. J. 648, 797 (2006), 

arXiv:astro-ph/0604485. 
[66] E. Komatsu, K. M. Smith, J. Dunkley, C. L. Bennett, 

B. Gold, G. Hinshaw, N. Jarosik, D. Larson, M. R. Nolta, 

L. Page, et al., ApJS 192, 18 (2011), 1001.4538. 
[67] S. Nesseris and L. Perivolaropoulos, JCAP 1, 18 (2007), 

arXiv;astro-ph/0610092. 
[68] L. Perivolaropoulos, Phys. Rev. D 71, 063503 (2005), 

URL http : //link . aps . org/doi/10 . 1103/PliysRevD . , 
I 71.0635031 

[69] J. M. Cline, S. Jeon, and G. D. Moore, Phys. Rev. D 70, 

043543 (2004), arXiv;hep-ph/0311312. 
[70] R. Gannouji, D. Polarski, A. Ranquet, and A. A. 

Starobinsky, JCAP 9, 16 (2006), arXiv:astro- 

ph/0606287. 

[71] M. C. B. Abdalla, S. Nojiri, and S. D. Odintsov, Classical 
and Quantum Gravity 22, L35 (2005), hep-th/0409177. 



[72] L. Amendola and S. Tsujikawa, Physics Letters B 660, 

125 (2008), 0705.0396. 
[73] J. Martin, C. Schimd, and J. -P. Uzan, Physical Review 

Letters 96, 061303 (2006), arXiv:astro-ph/0510208. 
[74] M. Libanov, V. Rubakov, E. Papantonopoulos, M. Sami, 

and S. Tsujikawa, JCAP 8, 10 (2007), 0704.1848. 
[75] M. Kunz and D. Sapone, Phys. Rev. D 74, 123503 

(2006), arXiv:astro-ph/0609040. 
[76] G.-B. Zhao, T. Giannantonio, L. Pogosian, A. Silvestri, 

D. J. Bacon, et al., Phys.Rev. D81, 103510 (2010), 

1003.0001. 

[77] S. Tsujikawa, Phys. Rev. D 76, 023514 (2007), 
0705.1032. 

[78] A. de Felice, S. Mukohyama, and S. Tsujikawa, Phys. 
Rev. D 82, 023524 (2010), 1006.0281. 



